Longitudinal relaxation and thermoactivation of quantum superparamagnets 
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The relaxation mechanisms of a quantum nanomagnet are discussed in the frame of linear response 
theory. We use a spin Hamiltonian with a uniaxial potential barrier plus a Zeeman term. The spin, 
having arbitrary S, is coupled to a bosonic environment. From the eigenstructure of the relaxation 
matrix, we identify two main mechanisms, namely, thermal activation over the barrier, with a time 
scale Aj -1 , and a faster dynamics inside the potential wells, with characteristic time A^ 1 . This allows 
to introduce a simple analytical formula for the response, which agrees well with the exact numerical 
results, and cover experiments even under moderate to strong fields in the superparamagnetic range. 
In passing, we generalize known classical results for a number of quantities (e.g., integral relaxation 
times, initial decay time, Kramers rate), results that are recovered in the limit S — * aa. 



PACS numbers: 03.65.Yz, 05.40.-a, 75.50.Xx, 75.50.Tt 
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I. INTRODUCTION 

Superparamagnets are nanoscale magnetic solids or 
clusters effectively described in terms of their net spin. 
They have an internal anisotropy barrier and several local 
equilibrium states. The spin is coupled to the environ- 
mental degrees of freedom of the host material (phonons, 
nuclear spins, electrons, . . . ). This coupling provokes dis- 
turbances in the spin dynamics producing jumps over the 
barrier, thus magnetization reversals, giving rise to the 
phenomenon of superparamagnetism. The prefix super 
stands for the typically large net spins of these systems 
(S ~ lO^lO 4 ). 

Single-domain magnetic particles are an example of 
such nanoscale solids, having a magnetic moment of a 
few thousand Bohr magnetons [l] . Due to their enormous 
spin these particles can be described classically @,03- An- 
other example is provided by single- molecule magnets, 
such as Mni2, Fes, or Nii2 Their net spin is S ~ 10, 
so that the classical description is no longer valid, and 
a quantum treatment is required. Still, these molecules 
comprise some hundreds of atoms, yielding an interplay 
of quantumness and mesoscopicity that has stimulated 
an active research in recent years. 

Classically, the spin surmounts the barrier AU by ther- 
mal activation. When AU/k-gT 3> 1 the characteristic 
time for the overbarrier process can be approximately 
written in the Arrhenius form t = A^ 1 cx ey^p{AU / k^T) . 
This is the rotational countenpart of the Kramers' rate 
in the theory of activated processes in translational sys- 
tems [H 0| . The classical dynamics of these systems has 
been studied extensively. In particular, the analysis of 
longitudinal response revealed that the one-mode relax- 
ation picture with a characteristic time A] -1 is not valid 
in general |ji|8|]. A satisfactory description, however, is 
provided by a two mode response, one mode correspond- 
ing to the overbarrier flux and the other, with a time 
scale A" 1 , describing the faster intrawell dynamics 

Quantum mechanically, the spin reversal can take place 
by thermal activation or tunneling (or a combination 



of both) m m ll'2l Il3j . Tunneling occurs from one 
side of the barrier to the other (Fig. between reso- 
nant states "coupled" by transverse fields (or high-order 
anisotropy terms). However in the superparamagnetic 
regime (T > 2 K) thermal activation controls the physics: 
out-of-resonance, the crossings are only driven by ther- 
mal activation, while tunneling in resonance has the effect 
of lowering the barrier a few states. 

The theory of quantum thermal activation has focused 
mainly on the zero, or very small bias limit |14lll5| . Then 
the relaxation becomes well described by one mode, the 
overbarrier process. The effect of bias fields has been 
only studied on the thermoactivation rate ^(|, but not 
on the susceptibility, beyond the weak-field regime. Thus 
studies of the full dynamical response including the ef- 
fect of external fields are demanded. This is the issue we 
address in this work: the relaxation of a spin with arbi- 
trary S in contact with a bosonic (phonon) environment. 
The dynamics is studied in the frame of quantum mas- 
ter (balance) equations. We concentrate on the case of 
longitudinal applied fields, giving a complete character- 
ization of the relaxation mechanisms. The discreteness 
of the levels plays an important role in some aspects of 
the dynamics. The connection with the classical theory 
is done taking the S — ► oo limit (whenever well defined) 
recovering known classical results. 

The paper is organized as follows. In Sec.HIIwe present 
the quantum balance equations. The tools of linear- 
response theory employed and the eigenstructure of the 
relaxation matrix are discussed in Sec. IIIII The relax- 
ation eigenvalue spectrum suggest the introduction of a 
two-mode description of the susceptibility (Sec. llVfl fol- 
lowing closely the approach of Kalmykov and co-workers 
in the classical limit 0, 0] . The comparison between ex- 
act numerical results and such analytical approximation 
is done in Sec. |VJ We also give a complete characteri- 
zation of the response and simplified expressions (both 
for the susceptibility and relaxation times) in the range 
AU/k-eT ~ 10-20, experimentally the most relevant for 
superparamagnets. Technical details for some calcula- 
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tions are sent to the appendices. 

II. QUANTUM SPIN IN A BOSONIC BATH 

A. Spin Hamiltonian 

The minimal Hamiltonian capturing the physics of su- 
perparamagnets includes a uniaxial anisotropy term plus 
the Zeeman coupling with external fields. Here we study 
the longitudinal relaxation. Thus only fields parallel to 
the anisotropy axis will be considered: 



n = -ds: - b z s z 



(i) 



In the standard basis |m) this Hamiltonian provides a 
spectrum e m with a double well structure with minima 
at m = ±S (for B z < D(2S-1); see Fig.EJ. The barrier 
heights are AU± = e TOb — e±s, where mp is the index 
corresponding to the maximum level [l8| . For fields 
B z > D(2S — 1) the barrier disappears. 
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FIG. 1: Energy levels for a spin S = 10 described by the 
Hamiltonian (0 with D = 0.5 at B z = 1. Considering that 
D is given in Kelvin, the value used is closed to Mni2. 



B. Balance equations 

In contact with a bath, the dynamical equation for the 
diagonal terms of the den sity matrix p typically has a 
balance-equation form EH |2£j : 



N m = R 



m\m + lNm+l 



P, 



[P, 



m\m- 



ra-J-1 |m-^m) 
P-m — 1 1 77i ) j 



(2) 



where N m = p mm — (m\p\m). The transition probabil- 
ities P m \ m i depend on the energy differences A mm i = 



m|m' 

e m i. They fulfill the detailed balance condition: 



P 



m\m' 



'P 



m' \mi 



(3) 



which ensures that the Gibbs distribution, Nm = 
e -m(m)^2 t would be a stationary solution of © (N m = 
0). The balance equations can be obtained in a phe- 
nomenological way: postulating the form and then 
calculating the probabilities P m \ m > by means of Fermi's 
golden rule [hH ]. 

Equation J2J is a set of linear differential equations. 
In matrix form it can be written as N = — 72. N, with 
[N] m = N m and 1Z the relaxation matrix. We assume TZ 
diagonalizable, and write the solution for N m as: 



2S 



d e 



-Ait 



|Ai) + |A ) 



(4) 



with Aj the eigenvalues of 1Z and |Aj) the associated 
eigenvectors. The last term in Q corresponds to the 
equilibrium solution, i.e., |Ao) m = Nm\ The relaxation 
matrix 1Z must give Re(Ai) > 0, ensuring that the asymp- 
totic solution will be the equilibrium one poj . 



C. Spin plus bath formulation 

Apart from phenomenologically, the balance equa- 
tion s |2l l can be obtained in a more rigorous way 
[ToL Il6l l20| . One starts from a total Hamiltonian rep- 
resenting the spin plus its environment: 

Wtat = H(S) + J2 VMS)(a+ + a- q ) + H B - (5) 



Here Hb — J2 q Uqa^a-q is a bosonic bath modelling the 
host material, the F q (S) the spin-dependent part of the 
interaction, and V q coupling constants. If the spin-bath 
coupling is weak, the equation of motion for the den- 
sity matrix p can be obtained within perturbation theory. 
In absence of transverse fields, the diagonal elements of 
the equation for p form a closed set of kinetic equations 
like @, with transition probabilities 



Pm\m' \P 



(0) 



The L m , m / are matrix elements of the spin-dependent 
part of the coupling F(S) and the rates W m \ m i = 
W(A mm i), are given in terms of the following universal 
function 19] 



W(A) = 



X A K 



1 



(7) 



Here we have considered that the spectral density of the 
bath, J{oS) — | J2 q l^q| 2 — w q)> ^ as *he form J{tu) oc 
uj k at low L). When k = 1 the bath is called Ohmic, while 
if k > 1 it is called super- Ohmic. 

Having molecular magnets in mind, we will use the 
following coupling HH El: F(S) oc {S z , S±}, then 

L m , m ±i = (2m±l)^/S(S + 1) - m(m ± 1) E3- In addi- 
tion the bath will be super-Ohmic (k = 3). This models 



3 



the interaction with 3D phonons considering one-phonon 
emission and absorption process, the relevant ones at low 
temperatures for single-molecule magnets. The modifica- 
tions required to include other structures of the coupling, 
or other baths, would only involve the P m \ m > and they 
are easy to carry out. 



D. Classical limit 

As we would like to make the connection of our results 
with known classical results, let us briefly consider the 
S — » oo limit of the balance equations @ ■ To this end it 
is useful to introduce the scaled quantities: 



PDS 2 



Z = PB Z S 



h = 



(8) 



In terms of them the scaled Hamiltonian reads (37i — 
—a(m/S) 2 — £(m/5). The first two quantities are equiv- 
alent to those used in the description of classical magnetic 
nanoparticles. The "reduced field" h, is B z measured in 
terms of the field for barrier disappearance D(2S — 1); it 
differs from the usual classical definition h c \ = £/ (2a), as 
a consequence of the discreteness of the energy levels. 

The quantities a and £ should be kept constant when 
taking the limit S — * oo. Then h — > h c \ and f3Tt — ► 
—az 1 — £z = u(z). Physically more and more lev- 
els are introduced, towards a continuum, while keep- 
ing the anisotropy barrier and Zeeman energy constant. 
In this limit the transition frequencies A m m ±i tend to 
zero. Thus, no relaxation is left in the classical case if 
P m \ m ± 1 (0) = (see App.0). Such is the case for a pure 
super-Ohmic environment. A well defined classical limit 
(■Pm|m±i(0) 7^ 0) is obtained considering an Ohmic bath. 
For instance, for coupling to electron-hole excitations or 
just adding two phonon processes (Raman scattering) in 
the interaction to phonons. In those cases, the balance 
equation © goes in the limit S — > oo over a partial dif- 
ferential (Fokker- Planck) equation: 



dW 
~~dt 



d_ 

dz 



D(z)( 



dW du 
V dz dz 



W 



(9) 



with W(z,t) the probability distribution of z (~ m/S) 
and P m+ i| m (0)/5 2 -> D(z) (see App.|X]for the details). 



III. ANALYSIS OF THE LONGITUDINAL 
RESPONSE 

The purpose of this section is to present the theoret- 
ical tools necessary to describe the relaxation of a spin 
S. To made a system to "relax" we must put it in a 
non-equilibrium situation (e.g., subjecting it to a pertur- 
bation in the external field). The response to this probe 
informs about the relaxation mechanisms of the system. 
We will use linear-response-theory tools [21j and the anal- 
ysis of the eigenstructure of the relaxation matrix. 



A. Linear response theory 

For the effect of the applied perturbation to reflect 
the intrinsic properties of the system the force should be 
suitable small. The weakness of the probe has the tech- 
nical advantage of allowing the use of linear response 
theory. Several "experiments" can be made to study 
the response, namely, subjecting it to a sudden constant 
"force" or by removing the force after having kept it on 
for a long time. One can also consider the response to a 
force oscillating with frequency f2. Linear response the- 
ory provides an intimate relationship between these types 
of measurements. 

Assume first that we have excited the system with a 
constant field 5B Z , keeping it on until the system is equi- 
librated (in a total field B z + SB Z ). Then, we switch the 
perturbation off and measure the response: 



AM z (t) = (M z (t)) - {M z ) , 



(10) 



with (M z )q the statistical average at equilibrium with 

B z = B° z and (M,(t)> = Yf m =-S mN ™{t)- Since the 
asymptotic solution for N is the equilibrium one in B z = 
B Q Z then (M z ) = (M z (t -» oo)). Next, introducing the 
solution Q for N m (t) in (M z (t)), we find 



2S 

AM z = X zSB z J2 a i e ~ 



Ait 



(11) 



where Xz is the equilibrium susceptibility: Xz = 
d(M z ) /dB z . In linear response Xz is equal to [(M z (0)) — 
{M z (oo))]/SB z . The coefficients ai are given in terms of 
the eigenvalues of the relaxation matrix 1Z and the coef- 
ficients Ci by [see Eq. (J2J]: 



1 3 
at = — } m\Ai) T 
Xz J-- 1 a 



(12) 



m=-S 



Here a has been redefined as c, = Ci/SB Z , to get rid of 
dependences on SB Z . The coefficients obey the nor- 
malization condition: Yh—i °i = 1 which follows from the 
definition ijT2|) and the equality: Xz = J2i J2 m m \^i)m- 
The Ci can be evaluated making t — in Eq. Q . The ini- 
tial condition is the system equilibrated at B = B Z +SB Z . 
The occupation numbers can be written in linear approx- 
imation as N m (t = 0) 
that the cVs obey 



a4 0) b1 0) 



8B z dB z N^ ] \ B o, 



so 



dm? 



dB z 



2S 

^ Ci | A^ ) m , 
i=l 



(13) 



which is an overdetermined set of 2S + 1 linear equations 
(there are 2S* Cj's). Both |Ai) and c\ can be obtained 
using standard numerical routines . 

Finally, let us consider the relation with the other re- 
laxation experiments. The case with the system equi- 
librated at B z — B z and a perturbation SB Z suddenly 



added is just the complementary situation to case dis- 
cussed above. More interesting is the oscillating field 
probe. Here we define the dynamical susceptibility, x(f2), 
as the coefficient which relates (in the frequency domain) 
response and excitation: AM z (f2) = %(f2)5I? z (f2) with 
g(£i) = f dt c~ lQt g(t). Then the response to a periodic 
perturbation 8B z (t) oc e~' ot is given by 

K ' V Jo AM z (0) J 

= *<£ttw <14) 

Equations and 1)14(1 provide the relation between the 
different relaxation experiments. The analysis of the re- 
sponse gives the time scales A" 1 and the weights dj's of 
the modes involved in the relaxation process. To con- 
clude, note that neither A, nor depend on the external 
probe but only on the intrinsic dynamics of the system, 
as it was demanded. 



B. Eigenvalue and eigenvector structure 

Attending at the evolution of N in Eq. we see that 
the eigenvalues of 1Z give the different time scales A^ 1 in 
the relaxation, while the eigenvectors provide the change 
in the population levels. The eigenvalues and eigenvec- 
tors for S — 10, obtained by numerical diagonalization 
of 1Z. are plotted in Figs. [3 and 
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FIG. 2: Eigenvalues of the relaxation matrix for S = 10. The 
spin-bath coupling is A = 1CF 9 , a — f3DS 2 = 15 and h — 
10" 3 . The arrow marks the lowest non- vanishing eigenvalue 
Ai. Inset: coefficients function of h for the same S, A, 

and a. 

Apart from the zero eigenvalue (not plotted in Fig.0), 
the eigenvalues correspond to two, well separated, sets 
of time scales. From the structure of the eigenvectors 
(Fig. |3J) we see that the dynamics induced by |Ai) pro- 
duces the increment in the population in one well and 
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FIG. 3: Eigenvectors of the relaxation matrix for S = 10, 
a = 15, h = 10~ 3 and A = 10~ 9 . At ft = the states A 2m -i 
and A2m (jn > 2) are degenerated. Numerical diagonalization 
gives even-odd eigenvectors, so we use a small h = 10~ 3 . 

the depopulation of the other. Since this transfer oc- 
curs across the barrier, |Ai) is associated to the overbar- 
rier dynamics. On the other hand, IA2} represents the 
transfer from m — to both wells or viceversa. Further, 
I A3) and IA4) account for dynamics involving m = ±1 
and m = ±S respectively, and so on. Eventually |A 17 )- 
|A 2 o=2s), involve the levels close the bottom of the wells. 
All these processes are called intrawell ones, since no 
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activation over the barrier is required. Thus the slow, 
and well separated, time scale Aj -1 corresponds to the 
overbarrier dynamics, whereas the set of (comparatively 
close) intrawell processes is responsible for the fast dy- 
namics. 

Let us briefly comment on two aspects of the eigen- 
structure. First, apart from Ao, Ai and A2, the remaining 
eigenvalues appear doubly degenerated (at ft = 0) [2j| . 
Physically each member of the degenerate pair describes 
identical processes in different wells, and at B z — both 
wells are equivalent. When B z ^ the equivalence is 
broken and the degeneration lifted. 

Finally a connection with the classical case can be 
sketched. The Fokker-Planck limit JHJ of the balance 
equations J2J) defines a Sturm-Liouville problem. From 
the corresponding theory we know that the eigenvalues of 
the differential operator are real and we can order them: 
AJ3 1 < A^ 1 < A^ 1 < . . . . Their corresponding eigenvectors 
|A^) have n nodes. Restricting ourselves to ft. = 0, the 
problem is invariant under the transformation z — > — z. 
Then there exists a basis where the states are even or 
odd under z — > —z. Then |Aq } (the Gibbs solution) 
is an even function with zero nodes. Next lA^ 1 ), would 
be an odd function since it has 1 node, lA^ 1 ) an even 
function (with 2 nodes), and so on. These features are 
retained in the discrete (5 < 00) case (Fig.|2Jl. Actually, 
from the degenerate states (e.g., IA3) and IA4)) we can 
form the symmetrical and antisymmetrical combinations 
I A 3 ) ± |A 4 ), fully recovering the even-odd picture. 



IV. ANALYTICAL APPROACH 

In general, the response will depend on 25 modes [see 
Eqs. ((TTj) and (fbEjl ]. However, taking into account the 
eigenvalue structure of the previous section, we approx- 
imate the response in terms of two main processes (see 
Sec. IIII 5|) . One accounting for the overbarrier dynamics, 
with characteristic time A^ 1 , and the other describing 
the intrawell dynamics with a time scale A" 1 . The latter 
is a kind of collective mode which describes the close set 
of intrawell processes. This motivation for the two-mode 
approximation is reinforced by the good results it yields 
in the classical limit 0, H3 • 



identified with known magnitudes of the problem. This 
can be carried out following the approach of Kalmykov 
and co-workers in the classical case @. They observed 
that we have two unknown parameters (a and A w ), while 
x(f2) can be evaluated in two extreme limits, namely, the 
low-f2 and high-f2 ranges. Then, taking these limits in the 
above two mode approximation, a and A w are identified. 
The low frequency behavior of x(^) is given by: 



(17) 



with Tj nt the integral relaxation time, defined as the area 
under the magnetization relaxation curve after a sudden 
change SB Z at t = 0: 



Tint 



dt 



AM z (t) 
AM z (0) 



(18) 



Next, in the high- 17 limit the susceptibility can be ex- 
panded as: 



x(n) = x. 



fir e f 



' c r 



d_ / AM z (t) 
dt \AM z (0) 



(19) 



t=o 



where r e f is the initial slope of the relaxation. 

The proposed formula l|16|) must obey the exact asymp- 
totic equations l|17fl and (|19fl . whence: 



Tint 
-el 1 



a/Ax + (1 - a)/A v 
aAi + (1 - a)A w , 



(20) 
(21) 



These equations can be solved for a and A w , yielding: 



TintT ef - 1 



A v 



AiTint - 2 + (AiTef)- 1 



AiTint - 1 



(22) 
(23) 



This reduces the problem to the obtainment of the three 
characteristic times Aj -1 , Ti nt , and r c f. Naturally, they 
can be expressed in terms of the relaxation eigenvalues: 



T int = > J -r 



ECLi 
. A~ 



= 2J a* A; 



(24) 



However, the goal is to bypass the eigenvalue computa- 
tion by calculating them directly. 



A. Bimodal approximation 

Suppose that the relaxation can be approximated by 
two exponentials: 

AM z (t) S X zSB z [«e- Alt + (1 - a) e"^ 4 ] , (15) 

or, in the frequency domain, 



Xz 1 + ifiAr 1 1 + iftAw 1 ' 

Note first that this bimodal description is exact for 5 = 1 
(see App. |P|) . In general a and A w are coefficients to be 



B. Calculation of Tint, Ai and r e f 

In the classical limit, Brown Q calculated the low- 
est eigenvalue Ai in the high barrier case (see also |24j). 
Besides Coffey and co-workers derived r e f (2f|. On the 
other hand, Garanin calculated Tint for a system of bal- 
ance equations to which Eq. (J2J can be reduced. In this 
subsection we are going to extend the results for Ai and 
T f to the quantum case (recovering the classical results 
in the limit 5 — > 00), and that of Tj nt to the generic sys- 
tem of balance equations. This will result in a closed 
expression for the bimodal formula i|16|) . 
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1. Integral relaxation time, r m t 

The calculation of Tj nt is based in that the susceptibil- 
can be calculated analytically up to first order in Q 
. Exploiting this fact, one finds for r int (see Add. |B*1: 



£ 

Xz 



4> 2 



s 



^ /V (0) P 



where 



(M z - j)N. 



(0) 



(25) 



(26) 



Notice that $ m only involves equilibrium averages, being 
independent of the spin-bath coupling model. 

The dependence of P m \ m i in the denominator on the 
energy differences, gives minima at the crossing fields |26| 
, which remind of the maxima in the relaxation rate due 
to tunneling. 



2. Lowest eigenvalue Ai 

The calculation of Ai, which corresponds to the 
Kramers rate, constitutes an important task by itself. 
Needless to say its relevance in the theory of activated 
processes [5, 6]. In the classical case it is possible to de- 
rive an expression for Ai 0, giving good results for not 
too strong fields |24|. In the quantum case, a closed an- 
alytical expression for Ai exists only for h = due to 
Villain and co-workers |14l and later on reexamined and 
improved by Wiirger [l5| . Their result can be written in 
a form appropiate for future comparison as: 



S-i 



' /V (0) P 

iv m -f>re+l 



(27) 



with 



= V 5 N {0) 



Both Ai from Eq. I|27(l and , together with the ex- 
act numerical Ai are displayed in Fig. 01 showing their 
agreement at h = 0. However, as we increment the field 
and Ai deviate from each other. This is natural, 
should t^I = Ai, the response would be described by 
one mode, namely the overbarrier [i.e., a = 1 in the bi- 
modal approximation, see Eq. 1201) ]. The same behavior 
was found in the classical limit pUg • The physical reason 
for the disagreement between Tj~ t and Ai is the intrawell 
processes entering into scene at finite h. Formally, at 
h <C 1, the quantity P m +i\ m Nm is highly peaked 
at the barrier (to = mb), and well approximated by the 
overbarrier process; then t^ 1 = Ai. Increasing the field, 
however, develops a second peak around the bottom 
of the lower well, to = S — 1, and the sum adquires a rel- 
evant contribution from the intrawell processes @ (see 
also Fig. |SJ|. 




-i -i -i 

T int ' 1 ' 



FIG. 4: Top: Comparison between the analytical r 5 ~ t 
[Eq. 1251 1. t^ 1 [Eq. 1271 1 and the numerical lowest eigenvalue 
Ai in the unbiased case as a function of a, the barrier height. 
The parameters are S = 10, T = 0.1 and A = 10~ 9 . Wiirger's 
result and ta have an maximum relative deviation of 10 -3 
in the plotted range (at the lowest barrier heights). Bottom: 
-- 1 [Eq. <E2J] and Ai for S = 10, A = 10~ 9 and 
15 as a function of the field. The relative difference be- 
tween and Ai never exceeds the 1% at not too strong 
fields (h < 0.5). For h > 0.7 we find numerical instabilities at 
these values of a and S. 



To generalise the calculation of Ai , we need an observ- 
able where only overbarrier processes are reflected. Then 
its associated integral relaxation time would probably be 
well approximated by Ai. This can be done generalizing 
to the quantum case the new integral relaxation time, ta, 
introduced by Garanin for the classical Fokker-Planck 
equation [27j ■ To this end, one chooses as observable the 
difference of the well populations A = N+ — N-. To 
obtain the integral relaxation time associated to A, we 
need its low-fi behavior, involving the "static" response 
XA ~ 9b z A [cf. the case of the magnetization fTTjl ]: 



Xa(«) 



dA(Q) 
dB z 



Xa(1 + i^T A 



(28) 
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This is done in App. [5] yielding the following ta 



V. ANALYSIS OF THE DYNAMIC RESPONSE 



TA 



3 iz} $ $^ 



with (rrib is the maximum level) 



In this section we study the full dynamical response. 

(29) 

We compare exact numerical results with the approxi- 
mate bimodal formula constructed in the previous sec- 
tion. We also present approximate tractable expressions. 



(0) 



[A - sgn(j - m b )] 



(30) 



It results that $ m <I>^ is nearly constant. Then the sum 
in Eq. (|29[1 has only one main contribution (around the 
barrier m = mb), not reflecting the dynamics inside the 
wells, as it was intended. 

We have seen in Fig. 0] that at h = 0, ta matches 
Wurger result (|27|l for Ai. In fact, it can be seen that 
our T/\\h=o reduces in the high barrier regime to Eq. (|27|l . 
(see App.jnj). Another limit where we have an expression 
for Ai is the classical one. Taking the limit S — > oo in 
ta we indeed obtain Brown's result in the high barrier 
range (App. Q 



(1-h) exp[-<r(l - h) 2 }}. 



h) exp[-o-(l + hf] 
(31) 



Finally, in the quantum case and h ^ 0, we have to check 
ta with the numerical results for A^ 1 ; this was done in 
Fig- El There one sees that, contrary to the ordinary inte- 
gral relaxation time, t^ 1 provides a remarkable approx- 
imation for the Kramers' rate. 



3. Effective time, r e f 

We conclude with the effective time. By definition T e f 
is the initial slope of the magnetization, see Eq. (|19|l . 
The initial condition is the system at thermal equilib- 
rium with B z — B z + 8B Z , so that N m (0) is known, 
whence N m (0) follows [directly from Eq. Next, 
dAM z /dt \t=o is calculated from N m (0) (see Add. IBIfor 
details). This gives the following expression for r e {' 







'ef 



m+1 m • 



(32) 



The classical limit of Eqs. (JUJ), and J^J give 

the expressions derived over the years from the Fokker- 
Planck equation JSJ). They are the sought analytical ex- 
pressions for the characteristic times Tj nt , Aj" 1 and T c f 
in the quantum case. Together with Eqs. 1221) and i|23|) 
they provide a closed formula for the response under the 
bimodal approximation, which can be checked against 
exact results. 



A. Exact results vs. bimodal formula 

To obtain numerically exact results for the response 
[see Eqs. (Illf) and (|14l) ]. we compute numerically the 
eigenvalues Aj and the amplitudes ai [see Eq. lfl3*)l]. Com- 
parison between them and the analytical formula i|16[l is 
shown in Figs. [5] and [fj] In Fig. [S] we plot the field de- 
pendence of the susceptibility spectra. The agreement is 
excellent in the low temperature range (large a). The 
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FIG. 5: Real parts (top) and imaginary parts (bottom) of 
the dynamical susceptibility of a spin S = 10 with A = 10 -9 , 
at T = 0.1 and a = 15 in various applied fields. The symbols 
represent the exact numerical results and the solid lines the 
bimodal formula 1161 . The dashed (thick) lines correspond to 
the bimodal formula but with the approximate a, Ai, and A w 
from Eqs. 1551. and lEol. 

(T-dependence is presented in Fig. |SJ The slow dynamics 
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10" 4 10" 2 10° 10 2 10 4 10 6 10 8 10 10 10 12 



QJA^ 

FIG. 6: Real parts (top) and imaginary parts (bottom) of 
the dynamical susceptibility of a spin S = 10 with A = 1CP 9 
at T = 0.1 and h = 0.2. Here a = 4, 6, 8, 10 and 12. 
Points represent the exact numerical results and lines are the 
theoretical formula 1161 . 

at Q ~ Ai is well described by the analytical expression, 
while reducing sufficiently a the bimodal approximation 
starts to fail at high frequencies £1 ~ A w . Specifically, 
the second peak in the exact imaginary part is no longer 
a single Lorentzian and a broadening with respect to the 
analytical curve is observed. 

The reason for the overall agreement between the bi- 
modal approximation and the exact results can be un- 
derstood as follows. The overbarrier dynamics must be 
well captured, since it is exactly given by Ai, while Ai is 
very well approximated by ta (see Sec. IIV|I . Next, in the 
two mode approximation we assumed that the intrawell 
modes are close to one another, so that they can be de- 
scribed with the collective scale A w . We have checked 
that the degree of closeness does not depend on er so as to 
explain the disagreement for low a and high in Fig. [5] 
What happens is that at high a the bimodal approxi- 
mation works very well because only intrawell processes 
involving the lowest well states are accessible, yielding 
a, ~ for the rest of the modes. For the example con- 



sidered, S = 10, one indeed sees that at a = 15, only an 
and a±s contribute (involving m = ±S and m ± (S — 1)), 
and to a smaller extent aig and dig (see the inset in 
Fig. Lowering enough the barrier more and more in- 
trawell modes play a role yielding a broadening of the 
curve for the imaginary part of x(f2). However, we have 
to go down to a = 6 to find large disagreements, while 
most experiments in these systems involve a > 10-20. 

B. Intrawell versus overbarrier processes 

Let us study the relative importance of intrawell and 
overbarrier dynamics in the response, a subject that re- 
cieved some attention in the classical case 0, 0. The 
competition between them becomes mediated by a in our 
Eq. (|16|) . In particular when a = 1 only relaxation across 
the barrier takes place while at a = the opposite oc- 
curs. From Fig-jSJandlBl we see that the relative weight of 
the two modes depends both on h and the barrier height. 
The field and a dependence of a is drawn in Fig. {7\ for 
different S. The main features are as follows. Concerning 
the h dependence, the transition from overbarrier dom- 
inated (a ~ 1) to intrawell dominated response (a ~ 0) 
is produced in a relatively small window of h. Secondly 
the transition occurs at larger values of h the smaller the 
value of S is. In addition, the a dependence of a shows 
that the value of S is also crucial in the limit a — > oo. In 
this limit a tends to 1 or to depending on S. 

Let us try to explain qualitatively these results, in- 
voking a statistical mechanical argument. Inspecting the 
level structure of the spin Hamiltonian, we see that at 
£ = <t/S the levels — m and m— 1 become degenerated. 
Such field corresponds to 

Then, if h < h* the first excited level is m = — S (in the 
other well), whereas if h > h* it is m = S — 1 (inside the 
same well). Approaching the large a limit we just need 
to consider the fundamental state and the first excited 
level. Then, when a — > oo, if h < h* intrawell process 
are inhibited, because only m — S and m — —S are 
populated, giving a = 1. On the contrary if ft, > h* , the 
states that contribute are m = —S and m = S — 1, in the 
same well, so that a = 0. 

The above argument also serves to explain the large a 
limits in Fig. [3 because the field used, h = 0.15, is lower 
than h*(S = 1,2,3), and greater than h*(S > 3). On the 
other hand, we have h* — > when S — > oo. This is in 
agreement with the classical calculation where a — > as 
soon as one makes h ^ in the limit a — > oo Even- 
tually our statistical mechanical argument also explains 
the S dependence of the onset of intrawell modes when 
increasing the field (Fig. EJ. Specifically, it gives a step 
function Q(h* — h) at a — > oo, while the remaining curves 
can be understood as a smoothing of the abrupt step due 
to finite temperature (Fig.0). 
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FIG. 7: Top: Field dependence of a± (points) and a (lines) 
for different values of S. We have fixed a = 15, A = 1CP 9 
and T = 0.1. Middle: Barrier dependence of m and a at 
the same S's. The effective field has been fixed at h — 0.15. 
Bottom: Field dependence of ai and the approximate a of 
Eq. (OHJ. The numb ers stand for the a value. Increasing a 
the curves approach a step function (solid line). At a — 30 
we find numerical instabilities so that for a = 50, 1000 only 
the approximate formula is plotted. 



C. Approximate expressions in the Kramers regime 

To conclude we present approximate formulas for the 
two time scales of the response, Ai and A w , in the high 
barrier (Kramers) regime, and for the parameter a con- 



trolling their relative weight in the response. In this 
regime a is well described by the formula (see Add. [Q: 



a = [1 + i (1 - qf e 



a(q+l) 2 (h-h*)-\ 



(34) 



where q = (S —1)/S. Returning for a moment to the pre- 
vious discussion, note that this formula captures the h, a 
and S dependences of Fig. {7\ for instance it approaches 
the step function Q(h — h*) when a — > oo, and gives its 
smoothening for finite a. Note finally that Eq. 134|l de- 
pends only on parameters of the spin Hamiltonian, being 
independent of the coupling (as in the classical case ||)- 

Let us go now into the Kramers rate. At low tem- 
peratures Ai can be expressed in terms of the barriers 
AU± = £ mb — £±s m a clear Arrhenius form (App.Q 



Ai = (e 



-/3AC/+ 



+ C 



-0AU- 



)( 



1 



1 



P, 



mb + l|mb 



P. 



mb — 1 1 mh 



Recall that the transition probabilities P ? 



(35) 



m\m' 

\L m .m'\ 2 W m \ m i [Eq. 10], involve the coupling matrix el- 
ements L m , m ±i = (2m ± l)^f S(S + 1) — m(m ± 1) and 
the rates W m \ m > — \ Af nm ,/(e l3Amm ' — 1). Equation (jSHJ) 
generalizes the zero-field result of Villain and co-workers 
|14j . In addition, with the appropriate L m>m ±i and the 
spectral index k (from J(w) — Aw K ), it is applicable to 
other couplings and baths. 

Finally we turn to the fastest mode A w . In order to fix 
some time scale in the problem, we measure A w relative 
to Ai. Under the same assumptions used to obtain the 
above a and Ai, we find for A w /Ai: 



Aw 
Ar 



1 



45 2 1 - a 



(36) 



Ps-i\s 



^ Pm h + l\m b 



+ 



1 



P_ 



S+1\S 



x 1 m b + l|m b 



rrib — 1 l^rib 

l 



P 



m b —l\m h 



While the parts involving the barriers AU± and a, are in- 
dependent of the coupling model, the transition rates at 
the wells and at the barrier are quite sensitive to the in- 
teraction. Thus the ratio A w /Ai can be used to compare 
different couplings and baths. 

Formulas (|34|l , l|35|) and (|36fl provide tractable expres- 
sions appropiate for the experiments in the superparam- 
agnetic range, out of resonance conditions (see App. El 
for further approximate expressions). Their reasonable 
agreement with the exact numerical results can be seen 
in Fig. [S] Under resonant conditions, we expect that our 
equations will still work after some modifications. For 
example, the barrier height will be lowered, since rrib 
would not be the barrier top but the state at which the 
activated tunnel process takes place. 
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VI. SUMMARY 

The relaxation theory of a quantum superparamagnet 
has been developed. We have studied the minimal model, 
namely a spin of arbitrary S with uniaxial anisotropy in 
a longitudinal field. Scarce results were available for the 
full dynamical response in the presence of external fields. 
Still, the popular single-Debye form for the susceptibil- 
ity spectra could be expected to fail, in analogy with 
the classical situation. Here, these topics have been ad- 
dressed in the frame of quantum balance equations rig- 
orously grounded on the system-plus-bath approach to 
quantum dissipative systems. 

We began analyzing the eigenstructure of the relax- 
ation matrix associated to the system of balance equa- 
tions. The form of the eigenvectors allowed us to iden- 
tify and classify the different relaxation mechanisms. In 
this way, full content has been given to popular state- 
ments about the relation of eigenvalues and relaxation 
modes. Besides, it has been put in connection with the 
Sturm-Liouville eigenstructure of the classical Fokker- 
Planck limit (parity, nodes, etc.). 

Two main processes emerge: activation over the 
anisotropy barrier, with a time scale , and a bunch 
of close (in logarithmic scale) fast intrawell modes, with 
an "average" time scale A" 1 . The identification and sep- 
aration of the modes suggest the introduction of an ap- 
proximate two-mode expression for the dynamical sus- 
ceptibility. Then, following the approach of Kalmykov 
et al. for the classical case, the parameters of that for- 
mula are expressed in terms of three characteristic quan- 
tities: integral relaxation time Ti nt , lowest non- vanishing 
eigenvalue Ai (Kramers rate), and the initial slope of the 
magnetization decay r e f. The usefulness of the classical 
approach was that these three times could be obtained 
analytically (which was done over the years for different 
relevant situations). Then, we were left with the task of 
finding them all in analytical form in the quantum do- 
main. This has been accomplished here, getting formulas 
for the quantum Tj nt , Ai (via the integral relaxation time 
for the population difference ta) , and r e f , which recover 
known classical results when taking the S — > oo limit. 
With them, one has a bimodal equation for the dynam- 
ical susceptibility of a spin with arbitrary S where all 
ingredients can be expressed in closed form. 

We have compared exact numerical results with such 
bimodal expression; the formula results to be quite 
accurate, specially in the superparamagnetic regime 
{AU/k-gT ~ 10-20; the range of major experimental in- 
terest). Its limits of validity have been assessed and in- 
terpreted in physical terms. Furthermore, in the range 
where the bimodal description works best, we have de- 
rived simple analytical expressions for Ai, A w and a (the 
parameter controlling the relative weight of the two effec- 
tive modes). These generalize formulas available at zero 
field, while they can also be applied to other structures 
of the coupling and spectral densities of the bath (e.g., 
Kondo coupling to electron- hole excitations). 



In the presence of terms not conmuting with S z in the 
Hamiltonian (e.g., transverse fields) the situation is ex- 
pected to be altered by tunnel events between resonant 
levels —m and m — k. Then, the balance equations loose 
the three-term recurrence form which allowed to derive 
the closed- form solutions [2^. However, in the super- 
paramagnetic range, thermal activation is still expected 
to govern the physics. Actually, in non-resonant condi- 
tions (the generic case), tunnel is inhibited, while at the 
resonances, it can be accounted for heuristically by lower- 
ing the effective barrier a few states. Thus, we hope that 
the two-mode picture will be of use with the appropriate 
amendments. 

The formulas derived cover experiments even under 
moderate and strong fields, where few work had been 
done. The onset of competence of the dynamics inside 
the wells with the overbarrier mode would occur at fields 
of a few Tesla, so that comparison with experimental data 
would be possible. We finally remark that the treatment 
employed covers from the deep quantum case (S ~ 1) to 
the classical regime (S ^> 1). Thus, the equations de- 
rived recover, in the limit S — > oo, the classical work on 
magnetic nanoparticles. 
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APPENDIX A: CLASSICAL LIMIT OF BALANCE 
EQUATIONS 

In this appendix we consider the classical limit of the 
balance equation J2J. Any limit procedure requires to 
specify which quantities are kept constant and which 
scaled variables are needed to monitor the evolution. 
For the classical limit, a natural choice is to maintain 
fixed a and £. At constant T this implies keeping the 
anisotropy-barrier height and amount of Zeeman energy 
constant. This means, that the levels tend to a con- 
tinuum, so A m , m= ti — ► 0. Then we need the behavior 
of P m \m±i near A m , m± i = 0. Expanding P m \ m > up to 
first order, with d x f(x)\ x= o = —d x f(—x)\ x =o, and the 
detailed balance condition @: 



dP, 



therefore, 



P Pm\m ! ( A m m > — 0), 



Pm\rn' — Pm\m'(fi) /?A m m r + C?(/3A m m , ) 



(Al) 
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Now, it results convenient to introduce the notation, where we have utilized (|A2|) . Now, writing p m -i = Pm + 

[Pm-i — Pm) and multiplying the right hand side of IJA3I1 
Pm = P m +i\ m (0); Pm-i = fm|m-i(0), (A2) by S 2 / S 2 we obtain, 



identifying through this appendix A„ 
inserting in Eq. @, 



/?A mm - and 



A„ 



JVtn-X 



(Prn-l + 



2 

Pm-l 



-i)iV„ 



(A3) 



S 2 



S 2 (AWi -2JV m + JV m _i) 



,1 



^S 2 S 2 J 



5 

S(A^ m - A m _i) + -A 

m,m — 1 



(A4) 



r 



In the first line of l|A4(l we recognize the usual discretiza- 
tion of the second derivative 



S 2 (N m+1 -2N m + N m - 1 ) 



d 2 W(z,t) 
dz 2 



(A5) 



where W(z, t) is the classical spin distribution. The third 
line can be identified with the discretizations of the first 
derivatives: 



SA 



m+l,m 



S(N m - N m - 
/3S (e m +i — e, 

Pm+l Pm 



s 2 



s 2 



dW(z,t) 

dz 
du(z) 

dz 
8D(z) 



dz 



(A6) 



here, we have utilized that p m /S 2 — > D(z). The 
function D(z) depends on the specific form of p m {= 
P m+1 | m (0)) and u(z) is the classical energy function 
[u(z) = j3TL(z)]. Finally, the second line is the discretiza- 

W(z, t), since 



s- 



1 



( 



-iJV n 



Using JUl, and jAT)l. we obtain the continuum 

version of ©, Eq. ©, where £>(z) = P m+ i| m (0) 

Let us emphasize that this derivation only makes use of 
detailed balance condition @ , and a nonvanishing transi- 
tion probability at zero frequency .P m +ii m (0), which has 



been tacitly assumed from the beginning. Considering 
interaction to electron- hole excitations (F oc S±), with 
an spectral density J(ui) oc u> then D(z) oc (1 — z 2 ) [see 
Eqs. © and (7J in Sec. \n\. On the other hand, cou- 
pling to phonons (including two phonon processes) yields 
D(z) oc z 2 (l -z 2 ). 



APPENDIX B: DETAILS FOR THE 
CALCULATION OF Tint, r A AND r cf 

Here we calculate the three time constants, Ti nt , t\ 
and r e f, which characterize the magnetization relaxation 
in the bimodal approximation [see Sec. II Vj . 

(%) Tint 

As Ti n t describes the low-f2 behavior of the susceptibil- 
ity our purpose here is to obtain x(f2) up to first order 
in £1 [sec Eq. i|17|)]. For future convenience, we write N m 
as, 



N m S N%> + N £) = N W (i + f3 qm ) 



(Bl) 



where N^} and N^} are the zero and first order (in 8B Z ) 
parts of the evolution for N m . Inserting (|B1|I in (J2J we 
obtain the set of equations necessary to calculate the re- 
sponse in linear approximation: 



P 



m-\-l\m 



P 



m—l\m 



here P„ 



P 



of R 



m+llm 



and P 



m|m' 



i^9m + -Pm+l|m(?m-|-l 

) (B2) 
= B2). To write (IB2H in terms 



m—l\m 



and the equality e 



the detailed balance equation © 
•Ac?, = A r m '' ) have been utilized. 
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We are interested in the low frequency behavior, so we 
expand 



(B3) 



We substitute l|B3jl in l|B2Jl . and we solve it perturba- 
tively in Q. The zero- frequency order gives: 



P , — P , — p , ( n y ' - n y { >) \ 

± m+l\m 1 m— l\m 1 m+l\m\Hm+l "rn ) 

1| 



+ Pm-Mmiq^-q^)- (B4) 



Since (|B4|I depends only on the differences <?„+i — Qm , 
and defining 



g = fl (°) _ o (0) 



(B5) 



equation (|B4fl can be casted in a two-term recurrence. 
With the "boundary" condition P_s-i\s — the solu- 
tion reads: 



permit to write q$ as 



^ = E 



m — l -r S j — 1 ~r 

<F, ^ ..(o) $ 



p ^(o) j p ^(o) 



(B13) 

where $i has been defined in (|26|) . The low-fi expan- 
sion (fT~7f) and IjBljl allows to express v mt as, 

iTint = — E mNX- (B14) 

Introducing (|B13|) in l!B14|l and using the equalities, 

s i-i , s-i , S 

2-^ 3 p ^(o) Z^ p ^(o) Z^ jv j 

j=-S i=-s i=-S ^i+l|i-' V i j=i+l 

(B15) 

and, 



9m = 1 

and accordingly to (|B5|) and (|B6|I . 

g<°> = m - M* 



(B6) 



(B7) 



where M z is the magnetization. To obtain 1|B7(1 the nor- 
malization condition: 



E N X = o. 



(B8) 



m— — S 



has been used. The above condition comes from 
J2 m = 0- Besides this follows from the unicity of 
the Fourier expansion of the norm and ^2 m N$h = 1 
plus J2m — 1 ■ Now, we write (|B2ll up to first order 
in Q, 

= ifigW +P m+ i\ m (q { ™+i -q£>) +^-1^(^-1 ~^)- 

(B9) 

Once more, <|B4|) depends only on the differences. Intro- 
ducing, 



' m — 1 m 



P m+1 \ m e-^(q^ +1 -q^), (B10) 



equation l)B9|) transforms onto a two term recurrence re- 
lation, obtaining for the r m elements, 



-m J2 q 

j=-s 



(Bll) 



where we have again invoked P_s-i\~s — 0- Now, (|B11|) 
with (IB 101) and the first-fi order normalization condition: 



-M z J2 Nf + £ J^ 0) = *, (B16) 

j=i+l J=i+1 

we finally obtain the formula l|25|l for the integral relax- 
ation time. 

(a) t a 

The time constant ta arises when instead of the mag- 
netization susceptibility, one considers the low-£! re- 
sponse of the difference of the well populations A 

s 

A = N+-N_= ]T sgn(m-m h )N m , (B17) 

m— — S 

with e mo the maximum level. Notice that the derivation 
of ta must follow the same steps that for the derivation 
of ri n t. Further taking into account the definition of A 
and from the expression for xa(^) in l|28|l ta is given in 
terms of qm (|B13f) as: 



ir A = — V sgn(m - m h )NWq$ (B18) 

m— — b 

Finally, using equalities l|B15|l and 

s s 
-A N f + E sgn(m-m b )ivj 0) =$f , (B19) 

j=i+l j=i+l 

where 3?^ is defined in l|30|) . formula (|29|l is readily ob- 
tained. 

(Hi) T c f 

First we notice that the definition for r e f (|19|) can be 
rewritten in terms of N m {t) as: 



E tfJW - o, 



(B12) 



7U — — S 



1 5 

— - > mN m (0). 



AM z (0) 



(B20) 



m=-S 
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N m (0), is obtained directly from the balance equation 
making t = in J2J. Now, we recall that the system is 
initially at equilibrium with field B z = B z + 5B Z . Up to 
first order in 8B Z this can be written: 



N m (0) = + 5B z j3N$(m - M z ). 
Next, we insert (|B21() in @. Considering: 

p , /v( 0) = p /v(°) 

r m\m±l 1 v m±l m±l|m iv m ? 

l|B20)l reads: 



(B21) 



(B22) 



I / y m|m— l^'m— 1 



X 



m=—S 



/ , m+l|m iv m I 



(B23) 



Finally, performing the change m — > m + 1 in the first 
summand we obtain the final formula 132fl for the inverse 
of the effective time. 



APPENDIX C: APPROXIMATE FORMULAS IN 
THE HIGH BARRIER CASE 

Here we derive approximate expressions for a, Ai and 
A w in the high barrier case. 

The technical dificulties deal mainly with the calcula- 
tion of the sums in Tint , ta and r c f . For that let us first 
consider the main approximation to handle with <£>^ and 
®m&m- I n the range considered <E> m can be well approx- 
imated by: 

$ m = $ B ; m < 5 - 2 

$ m = + m = 5-l (CI) 



where $ B S £™ b _ s (Af 2 ~ j) N j 0) = a,ui <]> 

E?=™ b +i( M * -i)ivj 0) = $ s -i - <&„ ib 

two main contributions in the sum for T m t, one around 
?7ib and the other one around m = 5 — 1 (see Fig. 
This last peak in the sum is due to $w- On the other 
hand the quantity $ m $^[ results well approximated for 
practical purposes by the constant: 



m b = 
Besides there are 



(C2) 



Then the sum 



with $ Ar = E™ b - s (A + W , -^ b - 

for ta has only one main contribution around mb (see 
Fig. EJ. 

Now, we give the main steps to calculate the formula 
(|34|l for a. In the bimodal approximation t ivA = a/A± + 
(1 — a) /A w . Then, and taking into account that A w 3> Ai 
and Ai = ta: 



Tint 
TA 



(C3) 




FIG. 8: Summands of T ln t 1251 1. — (0 ) — 



-, (solid line) and 



r A lUl, 5- 



! — m — (dashed line). Different h are drawn: 

h = 0.0.2, 0.4, 0.6. The rest of parameters are S = 10, a — 15, 
A = 1(T 9 . 

Next we use approximations (|C1|) and 1C2|I . Neglecting 
the second peak in the sum, which is produced by in- 
trawell dynamics (it would contribute to the summand 
(1 — a)/A w ) we write: 



a = 5 



XA 



(C4) 



here we have also used that in the high barrier limit 
$ B = S$ N (see Fig. |Hl . Notice that g2J depends only 
on equilibrium magnitudes. To compute Xz and xa we 
approximate the partition function by: 

Z= e'^coshf + e'^-^coshgf}, (C5) 

with q = (5—1) / S. This form for Z has been calculated 
considering only the states m = ±5 and m = ±(5 — 
1). Restricting ourselves to the lowest states, becomes 
justified when a 3> 1. However, at sufficiently high h, 
levels in the right well with m < S — 1 have less energy 
than m = —5+1 even m — —5, then they would be more 
populated. On the other hand, this is the minimal model 
which includes sufficient states to deal with intrawell and 
overbarrier contributions. Then, from l|C5|l we obtain \ z 
and xa- Finally, considering only the leading terms in a 
we obtain the final expression (|34J) for a. 

Now we turn our attention to Ai . In the main text we 
have argued that Ai = r^ 1 . Then, and considering l|C2|l . 
we write: 



TA 



XA 



S-l 

£ 

m=-S 



1 



N (0) p 
-i V rn 1 n 



(C6) 



+l|r 



We use now 8 m [see Eq. 
then it can be checked that: 

A £ 29 m 



i.e 



J2 S 



N (o) 

m+l iv j 



(C7) 
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and 



X a = 2d B J mh = 2/3$ J 



Thus A 1 1 reads: 



To obtain this formula we have used equation 1|C6() for Ai 
and (|C4|) for a together with $ B = S$ N . Now, taking 
(C8) only leading terms in the sum we readily obtain 11361) in 
the regime of interest. 



A r i -(i- 



s-i 



1 



iy(0) p 



(C9) 



m=-S JV '« ^m+ll 



First, Eq. HC9(I equals the Wiirger result l(T7Jl. since 
#m — $ro b = 1/2- Furthermore, the sums for 8 mb and 
l/A^ P m +i| m can be evaluated following Garanin in 
some limiting cases 16]. In particular the classical limit 
can be carried out. In this limit the sum is converted 
into an integral, which in turn can be solved. Then using 
approximate expressions for the classical Z when a 1 
(see Sec. II of [22j) and the same for 9 mb , we recover the 
Brown result (J3TJ for D(z) oc (l-z 2 ) [see Eq. @]. Even- 
tually, one also considers a/S > 1 to evaluate the sums 
in i|C9(l . In this case, becomes sufficient consider only the 
two main terms in the sum, namely 1 /N, 



(o) 



i/ivi 0) p ? 



Pm h + l\m b an d 

TOb _i| mb . In addition we only take into account 
the leading terms in a for mb , obtaining the Arrhcnius 
formula l|35|) for Ai. 

Finally we consider A w . In the expresion for A w 
[Eq.J2Hl)], the denominator can be approximated by 
A i Tint — 1 = a — 1 [see Eq. (IC3j) ]. Finally assuming 
t^ 1 > Ai (see Fig. EI) A w reads: 



APPENDIX D: EXACT ANALYTICAL 
SUSCEPTIBILITY FOR S=l 



When S — 1 it is possible to derive an exact analytical 
expression for x(^)- For that, we calculate the eigen- 
values of the, in this case, 3x3 relaxation matrix 1Z, 
obtaining: 



Ai,: 



(r + + r_) t \l (r+ - r_) 2 + ap^p%, (di) 



f±i|o[l 



plus the zero eigenvalue, Ao = 0. Here r± 
exp{/3(e±i — eo)}] Notice that the bimodal description 
is exact for S = 1 since A w = A2. Besides, using the 
relation of t^ 1 with the eigenvalues, i.e. t^ 1 — aAi + 
(1 - a)A 2 [Eq.||23] a yields: 



Then A w /Ai: 



A w ^ 



(CIO) 



1 



A2-A1 



(A 2 - r-/) , 



(D2) 



A v 



1 a 



Ai 45 2 1 - a 



s-i 



/ j Jv m r m+l\ 



m— — S 
/ S-l 



E 



" at(0) p , 



with r 



/3/( x ,Z)(P 1 ( " ) +P^ |0 ) 



,(0) 



The two modes A 1; A 2 and a are introduced in ltTH|) 
. (Cll) obtaining the exact expression for the response. 
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